Antiferromagnetic Excitations and Van Hove Singularities in YBa 2 Cu 3 06+i 
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We show that in quasi-two-dimensional d-wave superconductors Van Hove singularities close to 
the Fermi surface lead to novel magnetic quasi-particle excitations. We calculate the temperature 
and doping dependence of dynamical magnetic susceptibility for YBCO and show that the proposed 
excitations are in agreement with inelastic neutron scattering experiments. In addition, the values of 
the gap parameter and in-plane antiferromagnetic coupling are much smaller than usually believed. 

PACS numbers: 74.72.Bk, 75.50.Ee, 61.12.Bt 

The spin dynamics of the cuprate superconductors is 
still a subject of controversy. The spin-spin correlation 
function XPh{<l,(lz,w) = x' + ix", Q = (<?*><?»), usually 
probed by NMR, magnetic Raman scattering or, more 
directly, by inelastic neutron scattering (INS) experi- 
ments, is still not well understood. INS measurements on 



YBa 2 Cu306-fn;r(|YBCO) have been performed by a num- 
ber of groupsErQ and have lead to a consistent picture of 
the measured physical susceptibility XPh(<i, <?zj u) for q 
near the antiferromagnetic (AF) wave vector Q = (n,ir). 
It has been found that the dynamical AF spin fluctua- 
tions persist over the whole range of doping (0.4 < x < 
1) and the following remarkable features have been re- 
ported: 

(a) In the superconducting (SC) state xph(Q, fc, is 
limited by an energy gap, E g (Ref. |J). In heavily doped 
regime (0.65 < x < 0.92) E g is roughly proportional to 
the SC transition temperature T c , that is E g « 3.4 fcpT c . 
In the weakly-doped regime E g is very small and it is a 
much stronger function of T c (Ref. |]). 

(b) In addition, in highly doped YBCO the spin 
excitation spectrum is characterized by a sharp peak 
at some energy E r > E g , e.g., for optimally doped 
YBa 2 Cu 3 06.92, E r = 41 meV (Refs. For under- 
doped samples a broad resonance has been observed at 
20 < E r < 40 meV (Refs. ||,||). The resonance excita- 
tion is localized in both energy and wave vector and has 
Xpfe(Q) Izi E r ) oc s'm 2 (q z d/2) modulation, where d is the 
distance between two adjacent Cu02-planes in the unit 
cell of YBCO. 

(c) At ui w 50 meV there is a high energy cut-off, i.e., 
the intensity sharply decreases, regardless of the doping 
level. 

(d) In the normal state Xp/i(Qj q z , is characterized 
by a broad feature, spread practically over the whole en- 
ergy range, with a characteristic energy u> ~ 30 meV, 
approximately independent of the oxygen concentration. 

There has been a number of theoretical proposals for 
the explanation of features (a)-(d). It has been shown 
by Lu£3 that for d-wave SC state %"(q, q z ,w) should ex- 
hibit anomalous resonance peaks in q space. Feature (b) 
has been treated in terms of spin-flip excitations across 



the SC gapHii2 or as collective mode excitations in the 
particle-particle ch.annel.t3 

In this paper we show that the YBCO bilayer structure 
and 2D band topology, renormalized for AF interactions, 
with d-wave SC order parameter, leads to novel magnetic 
quasi-particle excitations which explain all of the features 
(a)-(d), observed in INS. 

We start by calculating Xpfc(Q> Qz> w ) for YBCO: wc 
consider a system of two 2D planes, corresponding to 
the CU2O4 bilayer in a unit cell of YBCO, with an ef- 
fective band dispersion: ^ = — 2t[cos(k x ) + cos(k y )] — 
At 1 cos(fc x ) cos(ky)— 2t"{cos(2k x )+cos(2k y )]±t±— fj,. Here 
index (i) corresponds to the electron bands with odd (o) 
and even (e) wave functions with respect to the mirror 
plane between the adjacent planes. We choose the pa- 
rameters t, t', t", t± and /1 by requiring that the shape 
and depth of the band and the size of the Fermi sur- 
face (FS) are approximately equal to those observed in 
ARPES measurements. According to Ref. the even 
band is nearly half filled regardless of the oxygen content, 
while at optimal doping the odd band is shifted upwards 
by 2t± w 100 meV and thus it is far away from half filling. 
This yields band parameters t = 115 meV, t'/t = —0.2, 
t"/t = 0.25, and \i = —228 meV (density of electrons 
rj(°) = 0.55 and n^> — 0.9). Although one usually in- 
terprets the parameters t, t' and t" in terms of tight- 
binding overlap integrals, here they should be regarded 
as results of our fits to the experimental data and as such, 
they should include some of the effects of mass renormal- 
ization. The odd band has a bifurcated saddle point at 
q = (vr,0), with = 30 meV and = 10 mcV 

(Ref. |l]); however, none of our conclusions are altered 
if one assumes an extended Van Hove singularity (VHS) 
instead. The FS of the odd and even bands are shown in 
Fig. § 

In BCS tiusacy the 2D Fermi liquid spin susceptibility 
is given byn3'E3 
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Af is the 



quasiparticle dispersion law in the SC state, f(E^ ) 
l/[exp(E®/T) + l] is thi 
ence factors are given by 



l/[exp(£j^ /T) + 1] is the Fermi function and the coher- 
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We assume that the SC gap parameter A^ has a simple 
[cos(fc x ) — cos(fcj,)]/2, although self- 



form A 



(i) 



consistent numerical solutions of the gap equation sug- 
gest that A k may have a somewhat different shape.Eil In 
Eq. ([l]) r is the phenomenological damping parameter. 
Accounting for the interaction between the susceptibility 
bubbles, we calculate-tha.spip susceptibility renormalized 
in the RPA manner :EjMBE|I 

X ^")(q, U ) = X < <J '>(q,a;)/[l + J (i '^q)x^'W)], (3) 

where J (e ' e )(q) = J\\ + Jj_, jO>°)(q) = J\\ - J±, 
j(e.o)(q) = J(°< e )( q ) = 0, J||(q) is the in-plane mag- 
netic interaction and Jj_(q) is the intcrplane magnetic 
coupling. We do not take into account the vertex and 
self-energy corrections which would further renormalize 
X ^")(q, W ). _ 

For the bilayer case the physical susceptibility, pro- 
portional to the INS measured intensity, is given by 

Em,n ex P{ ! *( z m ~ Zn)}x (m ' n} w ) (Ref- H)j where m 
and n are layer indexes. If coherence is preserved within, 
but not between the bilayers, this expression can be 
rewritten: 

XPh(<i,qz,u) = 

2[ X ( e ' e )(q, W )cos 2 (g z d/2) + X ( '°)(q, W )sin 2 ( gz d/2)]. (4) 

It is clear from Eq. (|4|) that the experimentally observed 
sm 2 (q z d/2) modulation of the INS (b) must originate 
from the odd band, provided the coupling between bi- 
layers is incoherent. 

For fixed q the quantity E^' 1 ^ (k, q) is a function of 
2D vector k and we denote its minima and saddle points 



J m%n ' 

and in the T, T 



by £^^(q) and £?|p (q) respectively. In a 2D system 



limit these extrema produce loga- 



rithmic singularities in either x'o (l> tJ ) or Xq~ % '^ (Qj w ) 
(see Eq. (Q)). For example, for p k " q 2 ^ 0, x'o (q, ">) 

diverges at w « I?^2(q), wm le Xo fa)'*') behaves like 
a step function. Alternatively, Xo(q> w ) diverges at the 
saddle point energies, Eip l \q), while Xo (q> w ) nas a 



"kink" .ta Provided that Jii and Ji arc not too large, the 
renormalized susceptibility (||) has poles (resonances) at 
u) close to both -E^V^(q) and (q). Physically, these 
poles in x^' 1 ' (Qj w ) describe spin-flip quasi-particle exci- 
tations (magnons). The self-energy corrections shift the 
position of the poles to the complex plane, producing the 
finite lifetime of the excitations, and more importantly, 
the logarithmic divergences become finite peaks, intro- 
ducing a minimal value of Jw) (q) for the creation of a 
quasi-particle. For quasi-2D system the resonances be- 
come broader due to the weak interaction between the 
bilayers. 

Low energy excitations at q = Q correspond to par- 
ticles and holes near the intercepts of the FS and the 
magnetic Brillouin zone (e.g., points A and B in Fig. 
0). The SC gap at these points A^ (= A^) is close 
to A^L, and thus the lowest energy excitation [feature 



(a) in INS] is E g = a£ + A£ = E™(Q). 
other hand, holes near VHSs and particles above the SC 
gap produce Eip'°\q) resonances. At q = Q this yields 
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the excitation energy E r — Esp'°\^; _i k - i k,-' 
where C and D are points depicted in Fig. |l| [feature (b)]. 
In Fig. U we present the energy dispersion of the quasi- 
particles, originating from both £^°^(q) and Esp°\q), 
for q along the high symmetry lines leading to Q, and 
assuming A^L = 15 meV. Away from Q the excita- 
tions clearly split into several bands, whose dispersion 
depends on the shape of VHSs, and £?|p (q) reaches 
approximately 0.3 eV near the magnetic Brillouin zone 
boundary. Of course, the excitations (resonances) shown 
in Fig. |2| may have very different lifetimes as well as in- 
tensities, depending not only on the FS geometry, but 
also on the appropriate coherence factors. 

Returning to the INS results, both E g and E r can be 
observed experimentally only if damping is sufficiently 
small. We calculated the single electron (hole) magnetic 
relaxation rate t^ 1 (see Eq. (9b) in Ref. ^6|) and found 
that in the SC state the relaxation is negligible for k 
near (w, 0) and energies below ~ 2.4A mox , in contrast 
with ordinary superconductors where Tj7 (u>, T) vanishes 
for lo < 3 A (Ref. |2^). As a result, any spin-flip excitation 
in x"^'^ (Q, ui), will not be damped for u smaller than 
even in d-wave superconductors. 



(o) 



E,° should be observed 



approximately 3.4A r 
Thus, both Eg and E r — A kc , ^ kc 
experimentally. On the other hand, in the normal state 
the E r excitation becomes heavily damped and is some- 
times referred to as the spin-pseudogap. We note that for 
the even band the saddle point energy lies far below the 
FS and the corresponding resonance excitation is over- 
damped even in the SC state and contributes only to a 
structureless background. 

The temperature dependence of the calculated odd 
band susceptibility x"^° ,0 \Q,^) is presented in Fig. ||a. 
We assume a BCS temperature dependence of the gap pa- 
rameter, although, admittedly, AtL(r) can be regarded 



2 



as a fitting parameter in order to match the experimen- 
tal INS spectraJfm Single band fits to NMRj£night shift 
suggest a larger value of SC gap parameterE3 than that 
we use here (AtL = 15 meV at T = 0). In the bilayer 
case only even band contributes to Xpfe(O) 0, 0) [Ref. p8[ , 
i.e., Aml, x may be larger than A^L (Ref. |||). In addi- 
tion, a different form of may alter Amax necessary 
to fit experiments. We assume J^°'°'(Q,) = —120 meV. 
The calculated spectra reproduce all of the observed INS 
features in the SC state: the gap originating from the SC 
gap parameter (a); the sharp resonance peak at 41 meV 
at low temperatures, due to spin-flip quasi-particle ex- 
citations related to the VHS in the band structure (b); 
the sharp drop of the signal above 45 meV, also due to 
the band topology and the effect of renormalization (c). 
The overall enhancement of X '°HQ> a; ) 1S due to the 

present in <i-wave 

At temperatures 



(o,o)2 



non- vanishing coherence factor p^ A £ 



superconductorsE (AJ^AJ^ < ) 
above T c the magnetic scattering is depressed. However, 
in agreement with the experimental data (d), a broad 
peak at w « 20 — 30 meV remains (spin-pseudogap) . Its 
origin is the same as that of E r resonance in the SC state 
(VHS); its position is shifted due to the SC gap vanishing. 

The doping dependence of ^"(°'°J (Q, cj) is shown in 
Fig. ||b. For the underdoped case we applied the same cri- 
teria as for the optimally doped, i.e., the band parameters 
are chosen in accordance with ARPES measurementstj 
which show that the VHS at 0) remains approxi- 
mately at the same energy below the FS, yet the filling 
is, obviously, much smaller. Therefore, in order to fit 
the experimental data we assumed, for the underdoped 
case, that t = 106 meV, t'/t = -0.18, t"/t = 0.18 and 
H = -125 meV (n^ = 0.7), with A^L = 5 meV and 
J(°>°)(Q) = -200 meV, and for the overdoped case we 
assumed t = 124 meV, t'/t = -0.22, t"/t = 0.24 and 

fj, = -204 meV (n (o) = 0.5) with A^L equal to that in 
the optimally doped case and J'°'°^(Q) = 0. The cal- 
culated result is qualitatively similar to the experimen- 
tal one.! Most importantly, since the value of A^L is 
smaller for the underdoped materials, but the distance 
from the VHS at q = (tt,0) to .the FS remains the same 
as in the optimally doped case,E3 the resonance E r must 
become broad at some doping level when this distance 
is approximately equal to 2AAmlix, and this is what the 
experiment suggests. 

In conclusion, we have shown that the Van Hove sin- 
gularities, observed in ARPES at approximately 30 meV 
below the Fermi level, lead to a novel, well-defined mag- 
netic quasi-particle excitations in the SC state, for quasi- 
2D systems and d-wave gap symmetry. Our model calcu- 
lations of Xp/i(Q> <7z, are consistent with all major fea- 
tures observed in INS experiments. The model suggests 
a relatively small value for A^L ~ l-9fcsT c , and a rela- 
tively weak Heisenberg interaction for the odd band of op- 
timally doped YBCO. In addition, our results show that 



the comparison of INS data and other magnetic probes, 
such as NMR, where both odd and even bands must be 
taken into account, may not be as straight-forward as 
previously believed. 

We are indebted to V. Barzykin, A. J. Leggctt, D. 
Pines and O. T. Vails for valuable discussions. This work 
has been supported in part by NSF cooperative agree- 
ment DMR 91-20000 through the STCS. 
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FIG. 1. Fermi surfaces of odd (solid line) and even bands 
(dashed line) in YBCO. The parameters are given in the 
text. The dash-dotted line shows the magnetic Brillouin zone. 
Points A and B (C and D) can be connected by the AF wave 
vector Q = (%, n) and correspond to E g (E r ) resonances ob- 
served in INS measurements (see text). 

FIG. 2. Magnetic quasi-particle dispersion for q along 
lines of high symmetry, in the vicinity of Q. The values are 
obtained from the minima (triangles) and saddle points (solid 
circles) of the particle-hole energy E+(k, q) [see Eq. (1)] at 
fixed momentum transfer q. 

FIG. 3. The odd band susceptibility x" ( °' o) (Q, w), denned 
in Eq. (3) and proportional to the INS measured intensity, as 
a function of energy. No finite lifetime effects are included 
(in both normal and SC state T = 0.5 meV). x" ( °' o) (Q,^) is 
gaped in the SC state and has a peak at Eip'°\ci) (see text) 
due to the VHS below the FS. (a) The results for optimally 
doped YBCO corresponding to T = 0, 80, 100, 200, 300 K 
(top to bottom), (b) The results for different doping levels 
and at T = K. 
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